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We propose a self-consistent model taking into account variations in adsorp- 
tion properties of the adsorbent surface in the process of adsorption-desorption 
of molecules of gas on it. We introduce a dimensionless coupling parameter that 
characterizes the interaction of an adsorbed molecule with polarized medium. It is 
established that the system can be bistable if the coupling parameter is greater than 
a critical value and the concentration of gas belongs to a certain interval. We show 
that the adsorption isotherms obtained within the framework of the proposed model 
essentially differ from the Langmuir isotherms and establish that the Zeldovich hys- 
teresis is possible. The kinetics of the surface coverage is analyzed in detail. We 
show that taking account of variations in adsorption properties of the surface in 
the process of adsorption-desorption leads to new phenomena: a "quasistationary" 
state in the case of the overdamped approximation and damped self-oscillations of 
the system in the general case. 

PACS numbers: 68.43.-h; 68.43.Mn; 68.43.Nr 



I. INTRODUCTION 

The study of adsorption of molecules on the surface of different bodies covers an ex- 
tremely wide class of problems of physics and chemistry and is one of the most important 
problems both from the theoretical point of view and for practical applications. The re- 
sults of numerous investigations show that adsorption of molecules on surfaces of bodies 
leads to changes in various physical and chemical characteristics of these bodies. The de- 



railed analysis o 



changes in the properties of the surface due to adsorption is given in 



LI, llfj- The results of investigation of changes in properties of the 
surface due to adsorption-desorption processes are also widely used in the design of various 
sensors (physical, chemical, and biological) [ll|, Q whose action is based on the use of the 
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change in a certain characteristic of a sensitive element of the sensor due to adsorption of 
molecules on its surface. 

The results of the theory of adsorption are extremely important for investigation of 
heterogenous-catalytic reactions because processes of adsorption and desorption are inte- 
gral stages of these reactions. 

The classical Langmuir theory that describes adsorption of a gas on solid surfaces is based 
on several assumptions. Numerous theoretical investigations, which, to a large extent, were 
stimulated by many experimental data that did not agree with conclusions of the Langmuir 
theory, were aimed at the construction of more general models free of one or several restric- 
tions of the Langmuir theory. An extensive material obtained on the basis of these models 
and applications to various problems of adsorption and catalysis are widely presented in the 
literature (see, e.g., 0, [l3, 14, Q, If], 0, [l8 ] ) . In particular, it is established that taking 
account of lateral interactions between adsorbed molecules can lead to a qualitative change 
in adsorption isotherms, namely, to a hysteresis of isotherms and to structural changes in 
ihe surface of bodies (surveys of theoretical and experimental results are given, e.g., in 

At the same time, as early as in 1938, in [l9j, Zeldovich has suggested an idea on a change 
in the surface in the course of adsorption and desorption due to the presence of adsorbed 
molecules on it. Using this idea, he has predicted a hysteresis of adsorption isotherms if 
the typical time of adsorption and desorption is much less than the relaxation time of the 
surface. 

Note that a change in adsorption isotherms due to lateral interactions between adsorbed 
molecules can also be interpreted as a consequence of a certain change in adsorption prop- 
erties of the surface caused by adsorption. However, as far as we know, the problem of 
variation in adsorption properties of the surface itself in the course of adsorption and des- 
orption of molecules of gas on it in the absence of interactions between adsorbed molecules 
and the possibility of hysteresis of adsorption isotherms in the this case remains open. 

The present paper is devoted to investigation of specific features of the behavior of adsorp- 
tion isotherms and the kinetics of the surface coverage by adsorbate molecules with regard 
for variations in adsorption properties of the surface in the course of adsorption-desorption. 

In Sec. 2, we propose a self-consistent model taking into account variations in adsorption 
properties of the surface in the process of adsorption-desorption of molecules of gas on it. 
We introduce a dimensionless coupling parameter that characterizes the interaction of an 
adsorbed molecule with polarized medium. We obtain adsorption isotherms and establish 
that their behavior essentially depends on the value of this parameter (Sec. 3). It is shown 
that, within the framework of the proposed model, the Zeldovich hysteresis is possible. In 
Sec. 4, we investigate specific features of the kinetics of the surface coverage. It is established 
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that variations in adsorption properties of the surface in the course of adsorption-desorption 
cardinally change the Langmuir kinetics. 



II. MODEL OF THE SURFACE WITH VARYING ADSORPTION PROPERTIES 

We consider a problem of adsorption of molecules of a one-component gas on the surface of 
a solid adsorbent. According to the classical Langmuir theory, molecules of gas are adsorbed 
on adsorption centers located on the adsorbent surface and the number of centers does not 
change with time. Furthermore, all centers have equal adsorption activity (energy-uniform 
surface), do not interact with each other, and each adsorption center can be bound only 
with one adsorbate molecule. The Langmuir kinetics of the quantity of adsorbed substance 
is described by the differential equation [l3| 

^ = k a c{i-e)-k d e, (i) 

where 9{t) = N(,(t)/N is the surface coverage by adsorbate, iV is the total number of 
adsorption centers, Nj>(t) and N (t) are, respectively, the numbers of occupied and free 
(N (t) = N — N b (t)) adsorption centers at the time t, k a and kd are the adsorption and 
desorption rate constants, respectively, and C is the concentration of molecules in the gas 
phase that is kept constant. 



The solution of this equation with zero initial condition #(0) = has the form 16j 



6(t) = e 

where 



t 

exp [ 

Tad 



(2) 



9? = -!— (3) 

is the stationary surface coverage (Langmuir isotherm), which is defined by the single di- 
mensionless quantity (dimensionless concentration) / = CK, K = k a /kd is the adsorption- 
desorption equilibrium constant for the given concentration of gas, 

Tad = KcTFd = TT1 (4) 

is the time taken for the surface coverage to reach the stationary value and = l/kd is 
the typical lifetime of a complex (adsorption center + adsorbed molecule). 

According to (2) and (3), there is a single- valued correspondence between the gas con- 



centration and the surface coverage. At the same time, in 19], Zeldovich has suggested an 



idea on a change in adsorption properties of the surface in the process of adsorption and 
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desorption and, using this idea, predicted a hysteresis of adsorption isotherms if the typical 
time of adsorption and desorption is much less than the relaxation time of the surface. 

In the present paper, within the framework of the Langmuir model, we take into account 
variations in properties of the surface of a homogenous adsorbent with plane boundary in 
the process of adsorption and desorption of molecules of a one-component gas on it. We 
introduce the Cartesian coordinate system centered at the surface of the adsorbent with 
OX-axis directed into the adsorbent perpendicularly to its surface so that the adsorbent 
and the gas occupy the regions x > and x < 0, respectively. Each adsorption center is 
simulated by a one-dimensional linear oscillator of mass m that oscillates perpendicularly 
to the surface about its equilibrium position x = 0. In the absence of adsorbate, the motion 
of an absorption center is described by the well-known equation of motion of a free linear 
oscillator 

m 0-77T + OL— + KX = 0, (5) 

at 1 at 

where k is the restoring force constant, a is the friction coefficient, and x is the coordinate 
of the oscillator. 

In the general case, due to occupation of the adsorption center with a molecule of adsor- 
bate, the electron structure of the center changes, which results in a change in the interaction 
of the center with neighboring atoms of the adsorbent, i.e., to a local polarization of the 
adsorbent. As a result, the medium acts on the complex with a certain force F p (f, t), where 
r is the running coordinate of the complex, which is the reaction of the medium on the 
electron reconstruction of the adsorption center. Under the action of this force, the complex 
tends to a new equilibrium position different from the equilibrium position x = of the free 
center. If the polarization of adsorbent caused by the formation of the complex is axially 
symmetric about the axis passing through the complex and parallel to the OX-axis, then this 
force has only the component normal to the boundary, F p (f, t) = e x F p (x, t), where e x is the 
unit vector along the OX-axis; for convenience, the center has the coordinate r— (x, 0,0). 

The force F p (f, t) either acts on each oscillator of the system if the oscillator is occupied 
with a molecule of gas or does not act if it is free, i.e., the oscillator interacts with the polar- 
ized medium only for discrete time intervals. Instead, we consider an approximation where 
the oscillator-medium interaction is continuous in time and the oscillator is permanently 
bound with molecule with the time-dependent probability equal to the surface coverage 9. 
In this case, F p (x,t) = F p (x) 9. This approximation is analogous to the mean-field approxi- 
mation used in problems of adsorption with regard for lateral interactions between adsorbed 
particles Q. 

We represent the force F p (x) in the form F p (x) = — dUi 2 x ^ , where Ui nt {x) is the potential 



5 



energy of interaction of the polarized medium with the complex. Expanding the quantity 
Ui n t{x) in the Taylor series and keeping only the linear term, we obtain 

U int (x) « -xx, (6) 



dUjntjx) 



is the force 

x=0 



where the parameter of complex-polarized medium interaction x — 
acting on the complex by the medium polarized by this complex. 

Ignoring the internal motion in the bound molecule-center system, i.e., considering the 
motion of the complex as a whole, and taking into account a change in the mass of the os- 
cillator in the process of adsorption-desorption within the framework of this approximation, 
we obtain the following equation of motion for the oscillator: 

where m e ff(9) = mo + m9 is the effective mass of the complex that varies in the process 
of adsorption-desorption and m is the mass of an adsorbate molecule. Since 9 < 1, the 
effective mass of the complex does not exceed its total mass M = mo + m = m e //(l). 

It follows from Eq. (7) that bonding of an adsorbate molecule with center leads to a shift 
of the equilibrium position of the oscillator by (%/ k) 9 and to a change in the potential energy 
of the free oscillator equal to kx 2 /2 by Ui n t(x) 9. In the limiting case where all centers are 
bound, i.e, 9 = 1, the equilibrium position x max = x/ K of the bound oscillator is maximally 
distant from the surface and the potential energy of the oscillator U(x) at this equilibrium 
position is minimal and equal to a half of the energy of interaction of the bound oscillator 
with the polarized medium, U min = U(x max ) = U int /2, U int = U int {x max ) = ~x 2 ''/ 'k. 

Thus, due to the interaction of adsorbate molecules with adsorption centers, centers shift 
relative to the nonperturbed surface of the adsorbent, i.e., this interaction leads to the 
formation (for adsorption) and healing (for desorption) of local defects of the surface. For 
Ui n t(x) < 0, these defects are "pits" (for x > 0) or "hills" (for x < 0) whose depth and height 
depend on the properties of both the adsorbate and the adsorbent. In the special case where 
all atoms of the surface are adsorption centers, this interaction leads to a shift of the surface 
of the adsorbent either inwards (for x > 0) or outwards (for x < 0), i.e., to the relaxation 
of the surface In other words, the processes of adsorption and desorption result in a 
deformation of the surface of the adsorbent, which leads, in the general case, to changes in 
the adsorption and desorption rates and, as a consequence, the surface coverage. Within 
the framework of the Langmuir theory of kinetics on the nondeformable surface (x = 0), the 
adsorption and desorption rate constants k a and ka do not depend on the concentration of 
gas and are defined by the Arrhenius relations 
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K = k + exp (- A.) , fc d = k_ exp (- A) , (8) 

where E a and ^ are the activation energies for adsorption and desorption, respectively, 
k + and are the preexponential factors, T is the absolute temperature, and ks is the 
Boltzmann constant. 

A molecule bound with center, due to its interaction with polarized medium, is in a deeper 
potential well than in the case of the nondeformable surface. Therefore, for its desorption, 
the molecule requires an energy greater than Ed by the value \ x i where x x = \Ui n t(x)\ is an 
additional energy that the bound molecule must acquire to break the bond with polarized 
medium. 

Generally speaking, the polarization of the medium can also affect the number of free 
molecules of gas that can overcome the adsorption barrier E a , i.e., a peculiar activation 
of free molecules of gas occurs and varies in the process of adsorption-desorption. Here, 
we do not take into account a change in the activation energy for adsorption (some results 
obtained with regard for a decrease in the activation energy for adsorption in the process 
of adsorption-desorption are presented in Appendix A). Supposing that the preexponential 
factor k^ is not changed, we obtain the following expression for fc d : 

k d {x) = k d esxp(--^^j. (9) 

It is worth noting that this quantity already depends on the concentration of gas because it is 
defined by the current state of the surface (the quantity x) that depends on the concentration 
of gas. Therefore, adsorption and desorption of molecules proceed on the surface whose 
adsorption characteristics vary with time. 

Introducing the dimensionless coordinate of oscillator £ = x/x max , we obtain the following 
autonomous system of nonlinear differential equations, which describes the kinetics of the 
quantity of adsorbed substance with regard for variations in adsorption properties of the 
surface in the process of adsorption-desorption: 



— = k a C(l-0)-k d 0exp(-g(), (10) 

where the dimensionless parameter g = \Ui nt \/kBT, which characterizes the interaction of 
an adsorbed molecule with polarized medium, can be called a coupling parameter. In the 
absence of interaction (the linear case) where x — (U int = 0), the parameter g = 0. 
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Note that system (10)-(11), in many respects, is analogous to the system of equations 
given in ^2(^ , which describes a transport of electrons in a system of molecules of biological 
nature with regard for electron-conformation interaction. 



III. STATIONARY CASE 

In the stationary case, it follows from Eq. (11) that £ = 9. Therefore, the equilibrium 
state of system (10)— (11) is defined no t by the pair of quantities (9 st ,£, st ), as is typical of 
dynamical systems of two equations 2l|, |22J, but only by one quantity 9 st , which is a solution 
of the equation 

I = F{6), where F{9) = — °— exp (-g0). (12) 

1 — 9 ^ ' 

In the general case, it is hardly possible to solve the transcendental equation (12) in the 
explicit form. Nevertheless, based on this equation, we can make a qualitative conclusion 
on the influence of a change in adsorption properties of the surface on the surface coverage. 
To this end, note that the ratio 9/(1 — 6) is equal to N^/Nq and the quantity I is also the 
ratio Nb/N but in the linear case. Rewriting relation (12) in the form 

^ = lex V (g9), (13) 



N 

we see that a change in properties of the surface caused by adsorption and desorption leads 
to an increase in the surface coverage for any concentration of gas. The difference between 
the numbers of bound centers in the nonlinear (g ^ 0) and linear cases increases with the 
coupling parameter g. 

This conclusion can also be made by taking into account that, in the stationary case, the 
desorption rate characteristic (9) has the form 

k d (9) = k d exp(-g9), (14) 

where the surface coverage 9 is a solution of Eq. (12). For a system whose adsorption prop- 
erties vary in the process of adsorption-desorption, the adsorption-desorption equilibrium 
constant 

K{9) = = K expM) (15) 

is greater than the classical adsorption-desorption equilibrium constant K. Therefore, the 
equilibrium of the system shifts towards an increase in the number of adsorbed molecules. 
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It follows from (14) that the interaction of adsorbed molecules with polarized medium 
results in an increase in the typical lifetime of complex 

r d (0) = = r d ex P (g 0) (16) 
kd{0) 

and, hence, an increase in the surface coverage. 

To analyze solutions of Eq. (12), we use a standard procedure 22j, namely: we consider 
the plane (0,1) and take into account that I > and < < 1 (Fig. 1). The required 
solutions of Eq. (12) are the abscissas of the points of intersection of a horizontal line 
corresponding to the given concentration I (the left-hand side of Eq. (12)) with the curve 
F(0) (the right-hand side of Eq. (12)), which is shown in the figure for different values of the 
parameter g. Since the behavior of the function F(0) is essentially different for g < g c and 
g > g c , where g c = 4, it is convenient to represent the parameter g in the form g = a g g c , 
where a g > 0. For g < g c (Fig. la), the function F(0) monotonically increases and lies 
to the right of curve 1 for the linear case (g = 0). Thus, as in the linear case, for any 
given concentration, the surface coverage has the unique value Of > Of, which agrees with 
conclusion made above on the basis of relations (13) and (15). With increase in a g , the 
curve F(0) becomes more deformed and its deviation from curve 1 increases, which leads to 
an increase in the difference Of — Of between the values of the surface coverage in nonlinear 
and linear cases. For a g = 1, the function F(0) (curve 1 in Fig. lb) has an inflection point 
for = C = 1/2 for the concentration I = l c = exp (—2) « 0.135. 

For a g > 1, the behavior of the function F(0) essentially changes: for the concentrations 
l\ and l\ (l\ < l\ < l c ) depending on the value of the parameter a g , the function F(0) has 
a minimum and a maximum at the points = 0\ > C and = 0\ < C , respectively, which 
are roots of the quadratic equation 

2 -0 + - = O (17) 
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and are equal to 




0\ 




The concentrations l\ and l\ corresponding to these surface coverages are defined as follows: 

l b n = (g0 b n -l)ex V (-g0 b n ), n=l,2. 

In Fig. lb, the concentrations l\ and l\ and the surface coverages 0\ and 0\ for them 
are shown with the use of dashed straight lines for the special case a g = 3/2. Thus, for 
concentrations l\ < I < l\, Eq. (12) has three solutions Of < Of < Of, furthermore, only 




FIG. 1: Graphical solution of Eq. (12) for different values of the parameter a g : (a) a g = (1), 0.1 
(2), 0.5 (3), 0.9 (4); (b) a g = 1 (1), 1.5 (2), 2 (3), 3 (4). Horizontal dashed straight lines correspond 
to constant values of the dimensionless concentration I. 
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FIG. 2: Bifurcation curve: branches 1 and 2 of the curve correspond to the bifurcation concentra- 
tions l\ and /rj, respectively. 



the first solution 9f lies near the linear With increase in a g , the concentrations l\ and 
l\ decrease and the difference between maximum and minimum solutions 9f — 9f increases. 

Analysis of the system of equations (10)-(11) shows that its stationary solutions 9f and 
#3* are asymptotically stable and the solution is unstable. 

If the concentration I tends to the end point of the interval (to the value l\ or Z^), 

then the stable 9f (or 9f) and unstable 9% solutions approach each other and, in the limit 
I = l\ (or I — I2), coalesce into one solution 9\ (or 6%) (in Fig. lb, for a g = 3/2, these cases 
are shown for curve 2.) Therefore, / = l\ and / = l\ are the bifurcation concentrations for 
which the dynamical system (10)— (11) is structurally unstable 2l|,|22| and has the compound 
(double) equilibrium states 9\ and 9\. These special cases should be investigated in their 
own rights. 

Using relations (12) and (17), we plot a bifurcation curve in the plane of parameters 
(a g , I). This curve defined in the parametric form 

a * = AO (1-6) ' 1 = Y=~9 eXp (-T^) (18) 

is shown in Fig. 2. For any point of this plane lying between the branches of the bifurcation 
curve, the system of equations (10)— (11) has three structurally stable equilibrium states: two 
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states are stable and one is unstable. If a point lies outside these branches, then the system 
has one structurally stable equilibrium state. At any point of the bifurcation curve, except 
for the cusp (a g — 1, I — l r ), the system has two equilibrium states: one is structurally 
stable and another is double [221 ] . At the cusp, the system of equations (10)— (11) has one 
triple equilibrium state [22 ]. 

The S'-shaped adsorption isotherm depicted in Fig. 3 for a g > 1 (curve 1) essentially 
differs from the Langmuir isotherm (curve 2) and, on the qualitative level, reproduces the 
Zeldovich hysteresis predicted in 19]. 

With quasistatic increase in the concentration from zero, the surface coverage, at the 
initial section of the lower stable branch OA of the isotherm, coincides with Langmuir one. 
For these concentrations, a released adsorption center manages to relax to the nonperturbed 
state before it binds with other molecule, furthermore, Td(0) ~ r^. With increase in the 
concentration up to the bifurcation value the difference between the typical lifetimes 
of complex Td{9) and T<f increases. In this case, a free adsorption center can bind with a 
subsequent molecule before it relaxes to the nonperturbed state. In this section of the lower 
stable branch of the isotherm, occupation of the surface by adsorbate molecules is determined 
by two factors: an increase in the concentration of gas and a change in adsorption properties 
of the surface. Due to the last factor, the isotherm deviates from the Langmuir isotherm, 
and this deviation increases with concentration. The pattern cardinally changes as soon as 
the concentration negligibly exceeds l\. In this case, the lower stable branch of the isotherm 
disappears and a new (unique) equilibrium state of the complex is considerably more distant 
from the surface than the previous one for I <l\. Furthermore, the passage to this state is 
performed for a constant concentration, i.e., solely due to a change in adsorption properties 
of the surface of adsorbent (according to the terminology used in 19], a slow adsorption 
occurs). This passage can require many molecules that successively take part in the process 
of adsorption-desorption on the same adsorption center. Thus, in this stage, a certain 
interaction between the molecule leaving the adsorption center and the molecule binding 
with it occurs. In Fig. 3, this stage of a sharp increase in the surface coverage for a constant 
concentration is shown by the dashed straight line AB. 

In passing to a stable equilibrium state lying on the upper stable branch of the isotherm 
(the point B in Fig. 3), the majority of adsorption centers is bound. As a result, a subsequent 
increase in the concentration of gas slightly affects an increase in the surface coverage. Such 
a "saturation" of the surface with adsorbate, which rapidly increases with parameter a g 
(Fig. 3b), occurs for concentrations considerably less than those in the linear case. 

In passing through the bifurcation concentration l\, the conditions of desorption for ad- 
sorbed molecules become essentially worse due to a considerable increase in the depth of 
the potential well and the displacement of the equilibrium state of bound adsorption centers 
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from the surface. As a result, for returning the system to the lower branch of the isotherm, 
the concentration should be considerably less than l\. With quasistatic decrease in the con- 
centration, the surface coverage decreases slightly and only, in approaching the bifurcation 
value l\, a variation in 9 becomes noticeable. In passing through the bifurcation concen- 
tration l\, the upper stable branch of the isotherm disappears and an equilibrium state of 
the complex lies considerably closer to the surface than the previous state for I > l\. As 
a result, the surface coverage sharply decreases for the fixed concentration due to a change 
(restoration) in properties of the adsorbent surface. The transition of the system from the 
upper branch of the isotherm to its lower branch is shown by the dashed straight line CD 
in Fig. 3. Note that this stage of drop of the quantity 9 is absent in 19]. A subsequent 
decrease in the concentration is accompanied by a decrease in the surface coverage along the 
lower branch of the isotherm, which, in fact, coincides with Langmuir isotherm. 

This behavior of the adsorption isotherm corresponds to the principle of perfect delay 23 ] 
according to which a system, which is in a stable state at the initial time, with variation in 
a parameter (concentration in the case at hand), remains in this state until the state exists. 

As the parameter a g increases, the bifurcation concentration l\ rapidly vanishes (see 
Fig. 2). Using the results of calculation, we can say that, for a g > 3, a change in adsorp- 
tion properties of the adsorbent in the process of adsorption-desorption leads to a peculiar 
adaptation of the system to a state in which the majority of adsorption centers are bound 
up to very low concentrations. 

Note that the isotherms obtained above for the surface whose adsorption properties vary 
in the course of adsorption-desorption (Fig. 3) are similar to the isotherms obtained with 
regard for lateral interactions between adsorbed molecules on a nondeformable surface 3, 
HQ and to the Hill-de Boer isotherms derived on the basis of the Hill-de Boer equation 
of state for adsorbed molecules (a two-dimensional analog of the Van der Waals equation) 



24]. 

In analysis of adsorption isotherms with regard for lateral interactions between adsorbed 
molecules (see, e.g., |25[]), for investigation of possible surface phase transitions, a critical 



temperature T c is introduced Bit 

I |24j . For the model considered in the present paper, 

using the expression for the coupling parameter g, the critical value g c , and the analysis 
of adsorption isotherms performed above, the critical temperature is defined as follows: 
ksT c = \Ui n t\/4:. For a system of adsorbed molecules, one stable state occurs for T > T c , 
whereas, for T < T c , two stable states are possible. The corresponding phase diagram for 
the adsorbed layer in the "surface coverage-critical temperature" coordinates is determined 
by the relation 
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L =A6 (l-9) 

and, as in the case taking into account lateral interactions between adsorbed molecules 
on a nondeformable surface within the framework of the mean- field approximation is 
symmetric about 9 = 1/2. 



IV. NONSTATIONARY CASE 



A. Overdamped Approximation 

First, we investigate the kinetics of system (10)-(11) within the framework of overdamped 
approximation where the masses of adsorption center and molecule are low and the friction 
coefficient is so large that the first term on the left-hand side of Eq. (11) can be neglected 
as compared with the second term. Using the well-known results for a linear free oscillator 
of constant mass 



2l| . this approximation is correct if 



r 2 M « r 2 r) (19) 

where tm = 1/lom, ^>m = \J k/M is the oscillation frequency of an oscillator of mass M, and 
r r = ol/k is the typical relaxation time of a massless oscillator. Since M is the maximally 
possible effective mass, condition (19) is even somewhat high. In this approximation, the 
system of equations (10)-(11) is simplified to the form 



- = l(l-9)-9exp(-gO, (20) 

where t' = t/ra is the dimensionless time and (3 = r T /rd ■ 

Analysis of system (20)-(21) performed on the basis of the qualitative theory of differential 



equations 22|, |26( shows that the stable equilibrium states of the system 9f and 9f are stable 
nodes and its unstable equilibrium state 9^ is a saddle. For the bifurcation concentration 
I = l\ (or 1 = 1%), the system is structurally unstable and has a compound equilibrium state 
9\ (or 9%), namely, a saddle-node with two saddle sectors and one stable nodal sector. The 
system is also structurally unstable for the critical concentration I = l c and a g = 1. In this 
case, the system has one equilibrium state 9 C , which is stable triple node. 

The numerical analysis of the system of equations with initial conditions for t' = 



9(0) = 0, f (0) = 



(22) 
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shows that, for any values of the parameters a g , /, and j3, the system monotonically evolutes 
to the nearest stable equilibrium state. Therefore, for the bistable system (a g > 1 and 
l\ < I < 1%), the stable equilibrium state 9f is inaccessible. The time taken for attaining the 
equilibrium state Of considerably depends on parameters, first of all, on the concentration. 

Let us investigate the kinetics of the surface coverage for a system that can be bistable 
for a g = 1.5. In this case, l\ « 0.0329 and l\ « 0.0754; 9\ « 0.789 and 9\ « 0.211. 

In Fig. 4, the kinetics of the surface coverage 6(t') is shown for concentrations less (/ = 
0.05, Fig. 4a) and greater (/ = 0.1, Fig. 4b) than the bifurcation concentration l\. For 
comparison, the Langmuir kinetics is shown in this figure by curve 2. For I < l\, the behavior 
of 9(f) is analogous to that in the Langmuir case: the quantity 8(t') monotonically increases 
from zero to the nearest stationary value 9f that lies near the linear value 9f (Fig. 4a). 
With increase in the concentration, this behavior remains true up to the bifurcation value 
l\ (moreover, both the stationary value 9f and the time taken for its attaining increase). 

For I > l\, the system has only one stable equilibrium state, furthermore, in this state, 
the surface coverage, which is close to the maximum possible value, is essentially greater 
than that in the linear case. Moreover, both the shape of the kinetic curve 9(f) and the 
time taken for attaining the stationary value considerably differ from the Langmuir ones 
(Fig. 4b). With increase in the concentration, the time taken for attaining the stationary 
value decreases. 

The value of the parameter f3 affects only the time taken for the system to attain the 
stationary value but does not qualitatively change the kinetics of 9(f) both for I < l\ and 
for I > l\. This time decreases if (3 decreases and increases if f3 increases, which is quite 
natural because a variation in f3 is equivalent to a variation in the friction coefficient a. 

For concentrations I > l\ near the bifurcation concentration 1%, the behavior of 9(f) 
qualitatively changes. The behavior of the quantity 9(f) for I — l\ (1 + 5), for low values 
of the relative concentration 5 = (l — 1%) /l\ > 0, is shown in Fig. 5. If the concentration 
slightly exceeds the bifurcation concentration l\ (curves 2 and 3), then the evolution of 9(f) 
can be conditionally divided into three stages: (i) from the initial zero value to a value of 
~ 9\ corresponding to the bifurcation concentration l\\ (ii) a very slow (in the limit lim , 

8— >0 

infinitely slow) variation in the neighborhood of 9\\ (iii) from ~ 9\ to the stationary value 9f. 
For low values of S, the time taken for attaining the stationary level 9f is determined mainly 
by the second ( "quasistationary" ) stage in which the system, in fact, does not change (curve 
2 in Fig. 5b). This behavior of the surface coverage 9(f) is caused by the well-known effect 
of slowing down of a system near a singular point for the bifurcation value of a parameter 



23j,|27j, 



2a . l29j in the case where a phase trajectory of the system moves near this point. 



This behavior of system (20)-(21) can be clearly explained on the basis of analysis of 
its phase trajectories in the phase plane (9, £). The phase trajectories of the system with 





FIG. 4: Kinetics of the surface coverage 6{t') on deformable (curve 1) and nondeformable (the 
Langmuir kinetics, curve 2) surfaces for a g = 1.5, f3 = 1; I = 0.05 (a), 0.1 (b). 
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FIG. 5: Kinetics of the surface coverage 9 it') for concentrations close to the bifurcation concentra- 
tion l\: a g = 1.5, /3 = 1; 5 = (1), 10~ 4 (2), 1(T 3 (3), 1(T 2 (4), lO" 1 (5). 
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FIG. 6: Phase trajectories of system (20)-(22) for concentrations less (I = 0.065) (a), equal(Z = l\) 
(b), and slightly greater (1 = 1% (1 + 6), 5 = 0.1) (c) than the bifurcation concentration l\\ a g = 1.5, 
(3 = 1. Dashed lines stand for the main isoclines of the system. 

zero initial condition (22) are shown for concentrations less (Fig. 6a), equal (Fig. 6b), and 
slightly greater (Fig. 6c) than the bifurcation concentration l\. The dashed lines in these 
figures stand for the main isoclines of the system: the isocline of horizontal slopes £ = 9 and 
the isocline of vertical slopes £ = (l/g) \n(9/l (1 — 6)). The points of intersection of these 
isoclines are singular points of the system. 

For I < l h 2 (Fig. 6a), the singular points A and C are stable (stable nodes) and the 
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singular point B is unstable (saddle). The phase trajectory in Fig. 6a starting from the 
origin of coordinates and going to the nearest singular point A completely lies between the 
main isoclines. Moreover, the immediate analysis of system (20)-(21) shows that all phase 
trajectories of system (20)-(21) with initial values belonging to the domain bounded by the 
sections of the main isoclines before their intersection at the point A also completely lie 
between the main isoclines. A change in the parameter (3 does not qualitatively change the 
behavior of the phase trajectories and only shifts them to one of the main isoclines: for 
(3 « 1 and (3 » 1, the phase trajectories are closely pressed to the isoclines of horizontal 
and vertical slopes, respectively. 

For I = l\, the singular points A and B coalesce into one (compound) singular point 
D = D{6\, £2), which is a point of tangency of the main isoclines (Fig. 6b). In this case, 
the phase trajectory is analogous to that in the previous case. For I > Z|, the system has 
only one (stable) singular point C (Fig. 6c). If the relative concentration is low, then a gap 
between the main isoclines in the neighborhood of their point of tangency D for I = l\ is also 
small. Since a phase trajectory does not leave the domain bounded by the main isoclines, 
it goes through the gap and, in a neighborhood of the point D, its motion becomes slower. 
Furthermore, the less the relative concentration, the narrower the gap between the main 
isoclines and the closer the phase trajectory approaches the point D and, hence, the more 
its slowing down near the point. This behavior of th e sy stem corresponds to the effect of 



critical slowing down near a degenerate critical point 23l. 1271. |28I|. 

As a result, for low values of 5, the function 9(f) in Fig. 5b has the form of a double step 
(curves 2 and 3). The first plateau of the step corresponds to the quasistationary state d\ 
and the second corresponds to the stable state Of. 

For (3 « 1, the kinetics of the surface coverage, which is shown in Figs. 4 and 5 for 
;he intermediate case (3 = 1, can be analyzed in a standard way with the use of a potential 



23j, |27(. It follows from Eq. (21) that £ = 0. Substituting this relation into (20), we obtain 



the following equation for 0: 

d0_ dV(0;l,g) 
M ~ dO ' 
where the potential V(0; l,g) can be represented as the sum 



(23) 



V(O;l,g) = V L (0;l) + V ind (O;9), (24) 

where 

l + l 
2 

is the parabolic potential for the Langmuir kinetics (1) and 



V L (O;l) = '-^O{e-20t) (25) 



20 



V ind (9; g) = j 2 {l - (1 + g 9) exp (-g 9)} - j . (26) 

is the potential caused by the action of the polarized medium on the complex. 

Thus, the behavior of the quantity 9 is completely defined by the form of the potential 
V(9; I, g) as a function of 9. Analysis of the potential V(9; I, g) shows that its form essentially 
depends on values of the parameters / and g. For g < g c , the function V(9] I, g), like Vl(9; I) 
in the linear case, has one minimum for a certain 9f. With increase in I, the value 9f 
increases and the minimum of V(9f t ; I, g) decreases. The behavior of the potential V(6; I, g) 
essentially changes for g > g c , l\ < I < l\. In this case, the function V(9; l,g) has the form 
of a double well with local minima at 9 = 9f and 9 = 9% separated by a maximum at 
9 = #2*. As the concentration I varies from l\ to 1%, the positions of the extrema, the depths 
of the wells, and the barrier between them 5V2 t i(l,g) = Z, — V(9f;l,g) vary. For 

concentrations near l\, the second well with minimum at 9 — 9% is rather flat, essentially 
shallower as compared with the first well with minimum at 9 = 9f, V(9f\ I, g) > V{9f, I, g), 
and corresponds to a possible metastable state of the system. As the concentration increases, 
the second well becomes deeper and the barrier between the wells decreases. At a certain 
concentration, the depth of the second well becomes equal to the depth of the first one. For 
higher concentrations, the second well is deeper than the first, V(9^; l,g) < V{9f] l,g), i.e., 
the system is in a metastable state in the first well and in a stable state in the second well. 
With a further increase in the concentration, the modulus of the difference between the 
minima of the wells \V(9^;l,g) — V(6f;l,g)\ increases and the barrier tfV^i^g) decreases 
(the slope of the first well between 9f and 9f is close to zero) and disappears for / = l\. As 
soon as the concentration becomes greater than the bifurcation value l\, the first minimum 
disappears, furthermore, for concentrations near l\ (/ > l b 2 ), the potential V(9;l,g), in a 
neighborhood of 9%, has almost a zero slope, which leads to the well-known critical slowing 



down of the system 23|, |27|, |28( . The behavior of this gradient dynamical system corresponds 



to the well-known principle of perfect delay 23f| according to which a transition of the bistable 
system between two stable states of equilibrium is absent. 



B. Influence of the Masses of Center and Molecule on the Kinetics of the Surface 

Coverage 

Now we investigate the kinetics of the surface coverage taking into account the masses 
of adsorption center and molecule. The estimating condition (19) for the overdamped ap- 
proximation can be represented in the form Q 2 M -C 1, where Qm = ^m/^Im — t mI t t and 
7az = a/2M are, respectively, the Q-factor and the damping constant of a free oscillator of 
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mass M. We also introduce the dimensionless quantities \i = m/m and d = (to/t^) 2 , 



where r = 1/u>q and uj = Jk/ttlq is the oscillation frequency of a free oscillator of 



mass m . The quantity Qm can be represented in the form Qm = Qo a/1 + A*, where 
(Jo = cuo/270 = r /r r = y/d/P and 70 = a/2m are, respectively, the Q-factor and the 
damping constant of a free oscillator of mass rriQ. 

Analysis of the singular points Of, Of, and Of of the system of equations (10)— (11) 
carried out in the phase space gives the following: For very low values of the Q-factor Qm, 
Of, Of and Of are, as in the overdamped case, stable nodes and a saddle, respectively [30] . 
With increase in the Q-factor Qm, starting from certain values that depend on values of the 
parameters a g ,l,/3, and /i, Of and Of become stable node-focuses. In the interval (0%, Of), 
near Of there appears a domain such that if Of falls within this domain, then the singular 
point Of is a saddle-focus, otherwise, it is a saddle. With increase in the Q-factor, this 
domain rapidly grows and covers the entire interval {0\, 0\) so that the unstable singular 
point Of is always a saddle-focus. 

Below, we give the results of numerical analysis obtained for a g = 1.1, i.e., for a system 
with possible bistability, and (3 — 1. 

The curves shown in Fig. 7 illustrate the influence of the masses of adsorption center and 
molecule on the kinetics of the surface coverage 0(t') for the concentration I = 0.11 lying 
in the middle of the bistability interval [If where l\ ~ 0.1064 and l\ ~ 0.1154. Here, 
for comparison, we present the quantity 0(t') without regard for the masses of adsorption 
center and molecule (curve 1) and the stationary level Of for this concentration (dashed 
line). For low values of masses (curve 2 in Fig. 7a), as in the overdamped case, the number 
of molecules adsorbed on the surface monotonically increases with time and reaches its 
maximum value Of. With increase in the coefficient d, which is proportional to the mass of 
adsorption center, the behavior of the quantity 0(f) changes. The surface coverage reaches 
its stationary value Of only after several oscillations about it (Fig. 7b). The amplitude and 
the number of oscillations as well as the time taken for attaining the stationary value Of 
increase with d (curves 2 and 3 in Fig. 7b). This behavior is caused by the inertia of the 
oscillator that overshoots its equilibrium position and the deviation of the oscillator from 
the equilibrium position increases with its mass. Therefore, taking account of the masses 
of adsorption center and molecule changes only the character of attainment of the nearest 
equilibrium state by the system. 

The behavior of the system can essentially change if the concentration is near the end 
point of the bistability interval. This case is shown in Fig. 8 where the stationary states of 
the system (the stable states Of and Of and the unstable state Of) are shown by dashed 
lines. The curves 0(t') in Fig. 8a show that, for low values of d, the behavior of the system, to 
a large extent, is analogous to its behavior for low masses considered above (Fig. 7a). Due to 




FIG. 7: Influence of the masses of adsorption center and molecule on the kinetics of the surface 
coverage 9(f) for the concentration I = 0.11 lying in the middle of the bistability interval; a g = 1.1, 
P = 1; m = 0, m = (curve 1), (j, = 1; (a) d = 1 (2), 5 (3), 10 (4); (b) d = 50 (2), 100 (3). 
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FIG. 8: Influence of the masses of adsorption center and molecule on the kinetics of the surface 
coverage 9(f) for the concentration I = 0.115 lying near the end point of the bistability interval; 
a g = 1.1, p = 1; m = 0,m = (curve 1), n = 1; (a) d = 10 (2), 20 (3), 30 (4); (b) d = 33 (2), 50 
(3), 100 (4); (c) d = 4900 (2), 5200 (3), 5300 (4). 

inertia, the system penetrates into the domain (9f, Of), which is a domain of attraction of 



the attractor Of (domain I) 22|. As in the overdamped case, with time, the system attains 
its stationary level 9f. However, as soon as the mass of the complex reaches a value for which 
the kinetic energy of the complex is sufficient to overcome the "barrier" 8f, the behavior 
of the system qualitatively changes (cf. curves 2 and 3 in Fig. 8b). Having fallen into the 
domain 6^*), which is a domain of attraction of the attractor 9f (domain II), the system 
moves toward its second stable stationary state Of. Having reached this state, the system 
oscillates about it with decreasing amplitude. As d increases, the time of attainment of the 



24 



e 1.0 




FIG. 9: Influence of the ratio of the masses of molecule and adsorption center on the kinetics of 
the surface coverage 9(f); a g = 1.1, /3 = 1, I = 0.115; mo = 0,m = (curve 1); d = 5300, fi = 0.5 
(2), 1 (3). 

stationary state 9f decreases because the system, in fact, is not delayed in the neighborhood 
of the unstable state Off (curve 4 in Fig. 8b). Thus, unlike the overdamped case for which 
the stable equilibrium state 9f is not attainable, due to the masses of adsorption center and 
molecule, the system can be in this state rather than in the state 9f. 

A further increase in d is accompanied by an increase in the amplitude of oscillations of 
the system about the stationary level 9f and the number of oscillations (high-Q system). 
However, the system does not leave domain II (curve 2 in Fig. 8c). If the amplitude of 
oscillations exceeds [Of —Off) , then the system falls into domain I and its subsequent behavior 
can be different. The system can remain in this domain and, after a time, it attains the 
stable stationary level 9f (curve 3 in Fig. 8c) as in the overdamped case (curve 1) and in 
the case of low values of masses (curves 2-4 in Fig. 8a and curve 2 in Fig. 8b). 

For a somewhat greater value of d, the kinetic energy of the complex can be sufficient for 
the system to overcome the "barrier" 9% for the second time. As a result, the system again 
penetrates into domain II and attains the stationary level 9% (curve 4 in Fig. 8c). It is clear 
that, with increase in mass (the value of d), the system can again return into domain I, etc. 
Therefore, the finial stable state of the bistable system (9f or 9^) depends on the number of 




FIG. 10: Influence of the masses of adsorption center and molecule on the kinetics of the surface 
coverage 0(f) for the concentration close to the bifurcation concentration a g = 1.5, (5 = 1, 
5 = 10" 4 ; m = 0,m = (curve 1); p = 1, d= 1 (2), 10 (3), 20 (4), 30 (5). 



26 



intersections of the unstable state Of by the system in the process of its evolution with time. 
The value of the parameter /3, which is proportional to the friction coefficient, also affects 
overcoming of the "barrier" Of by the system. Decreasing this parameter, i.e., increasing 
the Q-factor of the system, it is possible, in principle, to realize a mode in which the system 
visits each of two stable stationary states many times. This behavior of the system under 
study qualitatively agrees with the well-known behavior of a Newton gradient system whose 
potential energy has two minima separated by a barrier in the case of a low value of the 
damping parameter j^ . 




Since the Q-factor of the system is determined both by the mass of adsorption center and 
by the mass of adsorbed molecule, these two characteristics affect (but different in rights) 
the possibility of the system to overcome the "barrier" Of. The curves 0(t') depicted in 
Fig. 9, which describe the kinetics of the surface coverage for systems with equal masses 
of adsorption centers but different masses of adsorbed molecules, visually illustrate this 
conclusion. The behaviors of systems that returned from domain II into domain I are 
different. The kinetic energy of the complex with more light-weight molecule (curve 2) is 
insufficient for the complex to overcome the "barrier" Of for the second time, and the system 
is stabilized at the stable level Of. At the same time, the complex with heavier molecule is 
able to overcome this "barrier" and the system returns into domain II and attains the stable 
level Of (curve 3). 

Note that the behavior of the kinetic curves shown in Figs. 7-9, on the qualitative level, 
agrees with conclusion on the kinetics of the surface coverage made in Appendix B with the 
use of the effective potential in the special case where the relaxation time of the quantity 
0(t) and the characteristic times of the quantity £(£) are essentially different. 

For concentrations I = l\ (1 + S) that slightly exceed the bifurcation value l\, taking 
account of the masses of adsorption center and molecule can also qualitatively change the 
kinetics of the surface coverage (Fig. 10). For equal masses of adsorption center and molecule 
(/i = 1) and for values d < 1 (i.e., for Q < 1), the influence of masses on the behavior of 
0(t') is insignificant (curves 1 and 2, in fact, coincide). With increase in the mass, the plot of 
the function 0(t') in the form of a double step typical of the overdamped case remains true 
(curve 3 in Fig. 10b). However, the residence time of the system in the "quasistationary" 
state somewhat decreases, i.e., the delay of the system in a neighborhood of the singular 
point 0\ corresponding to the bifurcation concentration l\ is less than in the overdamped 
case. As the mass increases, the residence time of the system in the "quasistationary" state 
rapidly decreases (curve 4 in Fig. 10b) and, starting from certain values of d, the system 
moves to the stable stationary state Of without delay near the "quasistationary " state 0\ 
(curve 5). Thus, the possibility of an intermediate "quasistationary" state for the system 
essentially depends on relations between the inertial and dissipative characteristics of the 
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system. 

V. CONCLUSIONS 

In the present paper, we have proposed the self-consistent model taking into account 
variations in adsorption properties of the adsorbent surface in the process of adsorption- 
desorption of molecules of gas on it. Within the framework of this model, we have intro- 
duced the dimensionless coupling parameter that characterizes the interaction of an adsorbed 
molecule with polarized medium. We have established that the system can be bistable if 
this parameter is greater than critical and the concentration of molecules in the gas phase 
belongs to a certain interval. We have investigated bifurcation concentrations for which 
stable states of the system appear and disappear. We have obtained adsorption isotherms 
that essentially differ from the classical Langmuir isotherms. It is established the possibility 
of the Zeldovich hysteresis within the framework of the proposed model. It is shown that 
variations in adsorption properties of the surface in the course of adsorption-desorption can 
lead to a peculiar adaptation of the system to the state in which the majority of adsorption 
centers is bound up to very low concentrations. 

The detailed analysis of the kinetics of the surface coverage established that taking ac- 
count of variations in adsorption properties of the surface in the course of adsorption- 
desorption leads to new phenomena: a "quasistationary" state for the overdamped approx- 
imation and damped self-oscillations of the system in the general case. 
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APPENDIX: A 

A change in the activation energy for adsorption E a caused by polarization of the medium 
in the process of adsorption-desorption depends on various factors connected both with 
adsorbent and adsorbate. As an example, we consider the case where, due to polarization of 
the medium, the activation energy for adsorption E a decreases by the quantity x x equal to 
the increment of the activation energy for desorption E d caused by polarization. Supposing 
that the preexponential factor k + is not changed, we obtain the following expression for the 
adsorption rate characteristic k a (x): 
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k a (x) = k a exp [j^j , (A.l) 

which, like the desorption rate characteristic kd{x) defined by relation (9), depends on the 
concentration of gas. As a result, we obtain a system of equations that describes the kinetics 
of the quantity of adsorbed substance and differs from system (10)— (1 1) derived above under 
the assumption that the activation energy for adsorption does not vary in the process of 
adsorption-desorption only by the replacement of Eq. (10) by the equation 

— = k a C(l - 6) exp (g 0-k d 9 exp {-g £). (A.2) 

The factor exp (g£) in the first term on the right-hand side of this equation takes into account 
a change in the activation energy for adsorption in the process of adsorption-desorption of 
molecules of gas. 

In the stationary case, we obtain the same equation (12) for determination of the quantity 
6 st but with function F{6) with changed parameter g 

F(6) = T ^e X p(-2ge). (A.3) 

In addition, we obtain 

k (0) 

k a {6) = k a expM), K{6) = -^-L = K exp(2g6). (A.4) 

kd\y) 

Comparing (A4) and (15), we obtain a natural result that a decrease in the activation 
barrier in the process of adsorption-desorption leads to a shift of the equilibrium of the 
system towards an increase in the number of adsorbed molecules. 

Since the behavior of function (A.3) is identical to the behavior of function (12) with 
replacement of the critical value of the coupling parameter g c = 4 by g c = 2, the results of 
analysis of adsorption isotherms carried out in the the third section of the present paper 
remains also true in the case at hand with replacement g c — > g c /2. By analogy, the critical 
temperature below which the system can be bistable is changed, T c — > 2T C . 

The kinetics of the surface coverage 9(t') is analogous to the kinetics of 9{t') established 
above without regard for a change in the activation energy for adsorption (Figs. 4, 5, and 
7-10). However, in the considered case, it is somewhat faster, which is quite natural because 
a decrease in the barrier E a favors a faster filling of the surface with molecules of gas. The 
specific features of the kinetics of 6(t') depicted in Figs. 8 and 9 for a bistable system also 
occur, furthermore, they are realized for lesser values of masses of adsorption centers and 
molecules. 
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APPENDIX: B 

Here, we investigate the behavior of the dynamical system (10)-(11) that describes the 
kinetics of adsorption of molecules on the surface whose adsorption properties vary in the 
process of adsorption-desorption in the special case where the relaxation time of the quantity 
9(t) is much less than the characteristic times of the quantity £(£), i.e., the variables 9 and £ 
are fast and slow, respectively. Performing the adiabatic elimination of the fast variable 9(f) 



271 ] . namely, setting d6/dt = in Eqs. (10) and (11), we obtain the following representation 



for the surface coverage versus the slow variable £: 

The dimensionless coordinate of oscillator £ (t) is determined as a solution of the nonlinear 
differential equation 

d 2 £ d£ dU(0 
that describes the motion of the oscillator with effective variable mass 

rneffiO =m L + m9 L l^hg|l , {B .3) 

l + exp(-g£) 

m L = m + m9 L , (BA) 

in the effective potential 

Note that the second term in relation (B.3) for the effective mass disappears in the absence of 
polarization of the adsorbate in the process of adsorption-desorption, i.e., lim m e //(£) = mi. 

Therefore, we reduced the problem of investigation of the kinetics of the surface coverage 
to the problem of study of the motion of an oscillator of variable mass in potential (B.5). 
Since the quantity £ is the dimensional coordinate of a bound adsorption center, in terms 
of the coordinate x of this center, the equation of motion for it has the form 

. . d 2 x dx dU(x) 

where 

1 — exp (— bx 



m eff {x) 



- mi + m S + exp(-&*)' {B - 7) 
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FIG. 11: Normed effective potential for different concentrations: g = 4.4, I = 0.1 (1), 0.108 (2), 
0.11 (3), 0.115 (4), 0.117 (5). 
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Note that the effective potential (B.8) is analogous to the potential derived in the adia- 
batic approximation in [3JJ where the structural regulation of functioning of a macromolec- 
ular in repeating cycles of reactions is investigated. 

Analysis of the potential U(£) shows that, for g > g c , l\ < I < l\, it has the form of a 
double well with local minima at £ = £f and £ = £3* separated by a maximum at £ = 
where = 9^, n = 1,2, 3, and 9^ are the stationary surface coverages investigated in Sec. 3 
that satisfy Eq. (12). For g < g c and any I as well as for g > g c and / < l\ or I > l\, the 
potential U(£) has one minimum. 

The curves presented in Fig. 11 for g = 4.4 > g c clearly illustrate the essential influence 
of the concentration on the form of the potential. For concentrations lying outside the 



31 



interval the potential has a single minimum (curves 1 (for I < l\ ks 0.1064) and 5 

(for I > l\ as 0.1154)), furthermore, the equilibrium position of the oscillator for I > l\ is 
considerably more distant from the nonperturbed position £ = than that for I <l\. Curves 
2-4 illustrate the double-well character of the potential for the concentrations l\ < I < l\ 
and deepening of its wells (especially, the second well) with increase in the concentration. 

In the case of the double- well potential £/(£), the motion of the oscillator described by 
Eq. (B.2), which was initially at rest at the point £ = 0, can be different depending on the 
contributions of the inertial and dissipative terms. For small masses and large values of the 
friction coefficient, the oscillator rolls down into the nearest well of the potential [/(£) and, 
finally, is stabilized in the steady state at the point £f corresponding to a minimum of the 
potential. With increase in mass and/or a decrease in the friction coefficient, the kinetic 
energy of the oscillator may be sufficient to overcome the potential barrier between the wells 
and the oscillator falls in the second well. Depending on the values of masses and the friction 
coefficient, the oscillator can both remain in the second well with subsequent stabilization 
at its minimum and return to the first well. For very small values of a, the oscillator can 
many times visit each well before stabilization in one of them. 

The surface coverage 6(t) has a similar behavior. Therefore, the kinetics of the surface 
coverage shown in Figs. 7-9 and obtained without additional assumptions on fast and slow 
variables, on the qualitative level, agrees with conclusions made above on the basis of the 
motion of an oscillator in a double-well potential. 
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